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1.0  INTRODUCTION 


The  possibility  of  obtaining  improvement  in  the  ride  quality  of  a 
vehicle  using  a  damping  rate  which  varies  according  to  the  terrain 
roughness  has  been  considered  in  literatures  such  as  [  1  ]  -  [  5  ] . 

Recently,  a  preliminary  feasibility  study  on  the  identification  of  terrain 
roughness  and  frequency  characteristics  was  carried  out  by  Daigle  [  5  ] 
using  sampled-data  analysis  and  digital  filtering  techniques.  In  this  study, 
a  mathematical  decision  scheme  for  characterizing  the  terrain  roughness  in 
terms  of  its  frequency  (wave  lengths)  contents  is  developed.  The 
development  assumes  that  the  terrain  elevation  can  be  sensed,  sampled  and 
digitized  . .  The  terrain  elevation  sampled-data  stream  is  then  passed 
through  frequency  selective  filters  which  separate  the  frequency  components 
of  the  data  into  different  adjacent  bands  or  channels  on  the  frequency 
spectrum.  The  FMS  value  from  each  channel  is  computed,  and  the  relative 
amplitude  of  the  RMS  values  is  used  to  indicate  certain  degrees  of  terrain 
roughness  present  in  each  of  the  channels. 

Simulations  of  the  above  sampled-data,  digital  filtering  and  decision 
scheme  were  made  on  the  Systems  Engineering  Laboratories  (SEL)  digital 
computer.  The  effectiveness  of  the  scheme  for  indicating  or  identifying 
the  terrain  roughness  was  strongly  supported  by  the  simulation  results. 

The  use  of  the  terrain  roughness  identification  scheme  is  proposed  [  5  ], 
among  other  techniques,  as  a  possible  means  of  incorporating  a  microprocessor- 
based  on-board  adaptive  suspension  control  unit  for  a  vehicle. 


To  determine  the  feasibility  of  an  actual  implementation  of  the 
microprocessor-based  on-board  system,  a  preliminary  investigation  into 
experimental  microprocessor-based  filters  is  suggested.  A  main  concern  of 
the  investigation  is  the  computational  speed  and  numerical 
accuracies  of  the  microprocessor  in  the  realization  of  high  order  digital 
filters.  As  a  rough  guideline,  it  is  noted  that  the  digital  frequency 
selective  filters  used  in  the  formulation  of  the  above  terrain  roughness 
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identification  scheme  consist  of  bandpass  aiid  highpass  filters  whose 
critical  frequencies  are  less  than  10  Hz. 
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2.0  OBJECTIVES 


The  objective  of  this  report  is  to  investigate  the  design  and  to  carry 
out  the  actual  implementation  of  microprocessor-based  frequency  selective 
filters  which  may  be  suitable  for  the  terrain  roughness  identification  purposes. 
In  the  design  phase  of  the  investigation,  a  systematic  procedure  for  designing 
digital  filters  is  developed.  The  procedure  is  based  on  bilinear  transformation 
technique  with  emphasized  consideration  on  the  compensation  of  frequency 
warping  and  on  the  choice  of  the  ratio  of  working  frequency  to  sampling 
frequency.  The  effectiveness  of  the  proposed  design  procedure  will  be 
demonstrated  by  several  examples.  It  is  remarked  that  the  potential  of  the 
procedure  may  be  enhanced  by  the  incorporation  of  computer-aided  digital 
filter  design  techniques. 

In  the  implementation  phase  of  the  investigation,  the  hardware 
and  software  for  a  microprocessor-based  digital  signal  processing  system  will 
be  developed.  Microprocessor  realization  of  digital  frequency  selective 
filters  will  be  demonstrated  by  the  implementation  of  digital  lowpass, 
highpass,  bandpass  and  bandstop  filters.  The  actual  experimental  frequency 
responses  of  the  microprocessor-based  digital  filters  will  be  recorded 
and  compared  to  their  theoretical  frequency  responses. 

The  organization  of  this  report  is  as  follows.  The  systematic  procedure 
for  the  design  of  digital  filters  is  developed  in  Section  5.1,  the 

hardware  and  software  for  the  microprocessor-based  signal  processing  system 
and  filters  is  described  in  Section  5.2  and  Appendix  B.  The  actual  experimental 
frequency  response  of  the  microprocessor-based  digital  filters  is  given  in 
Section  5. ,3.  Section  4  discusses  the  results  of  the  investigation  and  provides 
a  few  recommendations  for  the  direction  of  future  effort.  A  summary  on  the 
design  of  analog  Butterworth  frequency  selective  filters  is  given  in 
Appendix  A. 
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3.0  CONCLUSIONS 


The  design  procedure  developed  in  Section  5  provides  a  systematic 
technique  for  obtaining  a  digital  filter  from  a  corresponding  analog  filter 
using  bilinear  transformation.  The  technique  takes  into  consideration  the 

frequency  warping  and  the  ratio  of  working  to  sampling  frequencies.  The 
effectiveness  of  the  design  procedures  is  illustrated  by  Examples  1-5, 
where  the  design  specifications  are  satisfactorily  fulfilled. 

Based  on  the  experimental  results  and  performance  of  the  microprocessor- 
based  frequency  selective  filters  presented  in  Section  5,  the  following  may 
be  inferred: 

.  The  computational  speed  of  the  microprocessor-based  system  is  sufficiently 
fast  for  the  implementation  of  the  digital  frequency  selective  filters  with 
the  required  specifications.  As  noted  in  the  Introduction,  the 

critical  frequencies  of  the  digital  filters  required  in  the  terrain 
roughness  identification  schemes  are  less  than  10  Hz.  It  is  seen  in 
Table  2  that  the  critical  frequencies  of  the  experimental  microprocessor- 

JU 

based  filters  can  be  much  higher  than  the  required  specification. 

This  further  implies  that  there  is  room  in  the  processing  time  for  implem¬ 
enting  higher  order  filters. 

.  The  numerical  accuracies  of  the  microprocessor-based  system  using  12-bit 
word  length  data  is  adequate  for  the  implementation  of  the  filters.  This 
is  clearly  illustrated  by  comparing  the  theoretical  frequency  responses  of 
the  digital  filters  depicted  in  Figs.  6,  8,  10,  12  and  14  to  the  actual 
experimental  frequency  responses  of  the  microprocessor-based  filters 
depicted  by  Figs.  19,  20,  21,  22  and  23. 


It  is  reminded  that  the  microprocessor-based  filters  can  readily  be  tuned, 
by  adjusting  the  sampling  frequency  ws ,  so  that  the  critical  frequencies 
coincide  with  the  desired  specifications. 
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4 . 0  RECOMMENDATIONS 


The  successful  preliminary  investigation  into  the  microprocessor  realization 
of  the  digital  frequency  selective  filters  provides  a  favorable  possibility 

for  implementing  a  microprocessor-based  on-board  terrain  roughness  identification 
system  using  digital  filtering  techniques.  The  following  effort  in  line  with 

the  investigation  of  microprocessor-based  signal  processing  system  in  this 
report  may  be  pursued  in  the  future: 

.  Implementation  of  microprocessor-based  system  with  parallel  processing; 

.  Use  of  16-bit  microprocessors; 

.  Use  of  fast  arithmetic  and  support  chips; 

.  Computer  aided  design  package  for  the  design  procedures  developed  in 
Section  5 

.  Implementation  of  the  mathematical  decision  criterion  for  identifying 
the  terrain  roughness  and  frequency  content  as  suggested  in  [5]. 

Some  of  these  efforts  are  currently  underway. 
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5.Q  DISCUSSION 


5.1  DESIGN  OF  DIGITAL  FILTERS 

The  basic  steps  in  the  design  of  digital  filters  generally  involve: 

the  specification  on  the  general  characteristics  of  the  filters; 
the  approximation  and  design  consideration  in  attaining  the  specific¬ 
ation; 

the  realization  of  the  filters  using  finite  precision  arithmetic. 

Step  (i)  depends  mainly  on  the  application  of  the  filters  while  step 
(ii)  depends  on  the  design  approach  adopted  by  the  designer.  Step  (iii) 
takes  into  account  the  limitations  of  digital  devices,  such  as  the  finite 
word  length  in  a  digital  circuit  or  machine  and  the  finite  computational 
speed. 

The  objective  of  this  section  is  concerned  with  step  (i-i)  of  the  design; 
step  (iii)  will  be  considered  in  the  next  section.  In  particular,  the 
design  of  infinite  impulse  response  (HR)  digital  filters  (lowpass, 
highpass,  bandpass  and  bandstop  filters)  using  analog  filter  formulas  and 
bilinear  transformation  will  be  presented  in  detail  in  this  section. 

While  a  digital  filter  may  directly  be  designed 
using  pole-zero  placement  technique  in  the  z-plane,  a  more  traditional 
approach  is  to  transform  an  analog  filter,  based  on  the  poles  and  zeroes  in 
the  s-plane,  into  a  corresponding  digital  filter  satisfying  a  prescribed 
specification.  Some  of  the  reasons  for  the  latter  approach  include  : 

.  the  straightforward  convergence  of  frequency  specification  (in  terms  of 
Hz  or  rad-s  *")  for  an  analog  filter  into  the  frequency  specification 
(in  terms  of  radian  frequency,  angle  around  the  unit  circle,  or  ratio 
of  frequencies)  for  the  digital  filter,  once  the  sampling  rate  is  given; 

.  the  utilization  of  the  highly  developed  art  in  the  design  of  a  variety 
of  analog  filters  to  obtain  the  corresponding  digital  filters  (e.g. 
Butterworth,  Chebychev  or  elliptic  filters) 

.  the  closed-form  design  formulas  for  analog  filters  which  can  be  translated 


(i) 

(ii) 

(iii) 
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to  yield  closed-form  design  formulas  for  the  corresponding  digital 
filters.  The  closed-form  formulas  facilitate  simplicity  in  the 
realization  of  the  filters. 

There  are  many  techniques  for  transforming  or  converting  an  analog 
filter  into  a  corresponding  digital  filter.  One  such  technique  is 
the  bilinear  transformation  which  is  described  below. 
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5.1.1  DESIGNING  DIGITAL  FILTERS  FROM  ANALOG  FILTERS  USING  BILINEAR 

TRANSFORMATION  WITH  WORKING  TO  SAMPLING  FREQUENCY  RATIO  CONSIDERATION 


The  design  of  digital  filters  from  application  of  bilinear  transf ormation 
to  the  formulas  of  analog  filters  has  been  considered  in  literature  such 
as  [  6  ]-[  7  ].  Most  of  the  design  procedures  in  these  literature, 
however,  do  not  include  a  systematic  way  for  determining  the  gain(say,  r  ) 
in  the  bilinear  transformation.  As  will  be  seen  shortly,  the  transformation 
gain  t  is  closely  related  to  the  quality  in  the  zero-order-hold  reproduction 
of  a  processed  sampled  data  from  an  analog  signal  by  a  digital  filter.  To 
ensure  a  desirable  reproduction  quality  in  the  digital  filter  output,  it  is 
important  that  a  proper  transformation  gain  is  used  in  the  design. 


In  this  section,  we  present  a  systematic  approach  to  the  design  of  digital 
filters  from  analog  filters  using  the  bilinear  transformation  which  takes 
the  ratio  of  working  or  critical  frequencies  of  the  digital  filter  to  the 
sampling  frequency  into  the  design  consideration.  The  approach  provides  a 
straightforward  procedure  for  choosing  the  transformation  gain  t  and  for 
obtaining  the  desirable  output  reproduction  quality  in  the  digital  filter. 

The  procedure  is  also  well  suited  for  use  in  computer-aided  digital 
filter  designs. 

BILINEAR  TRANSFORMATION 

For  sampled-data  signals,  the  Laplace  transform  (s-transf orm)  can 
be  shown  to  be  related  to  the  z-transform  by 


sT 

z  =  e 
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Cl) 


where  T  is  the  sampling  period.  Using  Pade’s  approximation  [  8  ],  (1)  can 
be  approximated  by 


sT 

e 


1  +  sT/2 
1  -  sT/2 


(2) 
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In  general,  one  may  redefine  the  mapping  as 

^  1  +  st/2 

z  =  - 

1  -  st/2 


2  z  -  1 


T  Z  +  1 


where  t  is  the  transformation  gain.  Relationship  (3)  is  known  as 
BILINEAR  TRANSFORMATION. 

The  mapping  of  the  s-plane  into  the  z-plane  by  bilinear  transformation 
(s)  is  shown  in  Fig.  1,  which  can  be  constructed  using  the  following 
relationships. 


Define  (see  also  Fig.  1) 

a)  =  analog  frequency  (rad-s  , 

^  ^  1 
a),  =  corresponding  digital  frequency  (rad-s  ), 

^  A  1 

a)g  =  2tt/T  =  sampling  frequency  (rad-s  )  , 

A  Wd 

0  =  —  2it  =  co.T  =  radian  frequency  (rad), 
co  d 

s 


(a)  Using  (3a),  the  frequency  axis  of  the  s-plane  (the  imaginary  axis, 

i  0 

s  =  jco  )  is  mapped  into  that  of  the  z-plane  (the  unit  circle,  z  =  eJ  ) 


as  follows: 


1  +  st/2 
1  -  st/2 


1  +  jc^T/2 
1  -  jcoaT/2 

+  (  co  t/2)2]  ej  tan-1(  V/2) 
a 


[1  +  U  t/2)  Vtan"  <-V/2) 
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(4) 

(5a) 

(5b) 


(b) 


The  relationship  (5)  represents  the  frequency  warping  or  distortion 
of  bilinear  transformation.  The  characteristic  of  the  distortion  is 
depicted  in  Figs.  2  and  3. 

From  (3a) ,  the  real  axis  (s  =  a)  of  the  s-plane  is  mapped  into  the 
z-plane  as  the  magnitude  of 


1  +  a x/2 

z  =  - 

1  -  ot/2 


(6) 


where  a  is  real,  and  where  it  is  seen  that 


-l<z<lfor-°°<a<0 


1  £  z  <  «  for  0  <_  a  <  2/t 
2 

-  co  <  z  <  _1  for  —  <  a  <  00 

T 

It  is  clear  from  the  above  that  the  left  half  of  the  s-plane  is  mapped  into 
the  unit  disk  of  the  z-plane. 


Now,  let  G(s)  denote  the  transfer  function  of  an  analog  filter  and  G(z) 
denote  that  of  a  corresponding  digital  filter.  Then  using  bilinear  transformation 
(3),  the  digital  filter  can  be  obtained  as 
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G(z)  =  G(s)  .  (7) 

2  z  -  1 

S  T  Z  +  1 

By  the  mapping  of  the  bilinear  transformation  (Fig.  1),  all  the  stable  poles 
of  G(s)  will  be  converted  into  stable  poles  in  G(z) .  Consequently,  the 
bilinear  transformation  (7)  always  yield  stable  digital  filters  from  stable 
analog  filters. 

The  effect  of  frequency  warping  or  distortion  (5)  on  the  analog  to 
digital  filter  conversion  (7)  is  illustrated  in  Fig.  3.  The  figure  also 
clearly  reveals  an  explanation  for  the  phenomenon  of  aliasings. 

It  is  remarked  that  the  transformation  gain  t  for  the  bilinear 
transformation  (7)  has  not  been  specified.  A  systematic  technique  for 
determining  t  is  given  in  the  sequel.  The  technique  also  automatically 
compensates  for  the  frequency  warping  or  distortion. 


DESIGN  CONSIDERATION 

The  digital  zero-order-hold  (ZOH)  reproduction  of  an  analog  signal 
having  a  dominant  working  frequency  «  (or  correspondingly  a)  )  depends  on 
the  ratio  Fig.  4  illustrates  the  variation,  with  respect  to  the 

ratio  o)^/cjs,  by  a  sample  and  ZOH  scheme!  As  can  be  seen  from  the  figure, 
the  "quality"  of  the  digital  reproduction  of  the  analog  signal  improves 
with  lower  ratio  of  u)^/ws.  It  is,  therefore,  desirable  to  design  a  digital 
filter  whose  working  or  critical  frequencies  are  much  lower  than  the 
sampling  frequencies.  A  first  design  consideration  in  a  digital  filter 
design  is  to  ensure  that  the  ratio 

0  <  oj , /oj  <<  .5  .  (8) 

d  s 

The  microproessor-based  sample  and  ZOH  scheme  is  described  in  Section  3. 
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Remark  1;  For  simplicity  and  clarity,  the  above  argument  is  approached 
from  time-domain  point  of  view  using  visual  experimental  results.  It  may 
be  remarked  that  similar  conclusions  can  be  obtained  using  frequency  domain 
analysis  [  6  ].  One  also  notes  that  (8)  is  in  agreement  to  the  Sampling 
Theorem  due  to  Shannon  and  Nyquist  [  6  "] . 

It  is  also  important  to  observe  the  time  delays  in  the  digital  outputs 
in  reference  to  the  continuous  signals  in  ^ig.  4. 


A 


Once  the  ratio  has  been  selected,  one  may  define  a  factor  R  as 


tan 


to 

s 


0  <  to^/cog  <  .5 


CO  T 
a 

2 


(9) 


where  (9)  follows  from  (5).  Note  that  0  <  R  <  °°.  T^rom  (9),  the  transforma¬ 
tion  gain  x  is  obtained  as 


T  = 


2R 

to 

a 


(10) 


Using  (8)  in  (3b)  ,  the  bilinear  transformation  becomes 

“a  z  -  1 
S  "  R  z  +  1 


(11) 


The  above  design  consideration  of  first  specifying  a  desired  ratio 
ud/co  thus  leads  to  a  systematic  choice  of  the  transformation  gain  t  for 
the  bilinear  transformation  as  shown  in  (11). 


of  an 


Using  (11)  as  the  basis  for  bilinear  transformation,  the  conversion 
analog  filter  G(s)  with  working  or  critical  frequency  to  to  a 
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corresponding  digital  filter  G(z)  with  working  or  critical  frequency 
0) , ,  follows  from  (7)  as 


G(z) 


G(s) 


s 


a  ~  i 
a  z— 1 

R  z+1 


(12) 


The  bilinear  transformation  in  (12)  ensures  that  the  desired 

co,/o3  will  be  obtained, 
d  s 

Finally,  it  is  important  to  note  that  the  working  or  critical 
frequency  w,  of  the  digital  filter  can  be  varied  by  simply  adjusting  the 
sampling  frequency  ui  . 
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z  -  plane 


of  s  -  plane  into  z  -  plane  by  bilinear  transformation 


Frequency  warping  in  bilinear  transformation 
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Fig. 4. 


“d/^s  =  1/15 


a)d/UJs"'=  "1/2 


ud/us  =  l 

ZOH  Reproduction  of  Sinusoidal  Signal  with  respect  to  Ratio 
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5.1.2  DESIGN  OF  DIGITAL  LOW? ASS  FILTERS 


The  design  of  an  n-th  order  analog  Butterworth  lowpass  filter  G^p(s) 
with  critical  cut-off  frequency  to  is  given  in  Appendix  A2.  To  obtain 

3.C 

a  corresponding  digital  Butterworth  lowpass  filter  G^p(z),  the  following 
systematic  procedure  may  be  used: 


(a)  Select  the  ratio  (o^/ojg,  where  to^c  is  the  desired  critical  cut-off 
frequency  of  the  digital  lowpass  filter  (see  Fig. 3a). 

(b)  Obtain  the  factor  R  and  transformation  gain  as 


R  =  tan 


o 


£0  T 

ac 


2 


(13a) 


or 


(13b) 


(c) 


Using  the  substitution  described  by  (12) ,  a  corresponding  n-th 
order  digital  Butterworth  lowpass  filter  is  obtained  as'*' 


glp(z) 


(s) 


S' 


to 

ac  z  -  1 
R  z  +  1 


(z+l)n 

=  K  - 

(z-p1)(z-p2)  ...  (z-Pn) 


with 


=  KG(z)  , 

A  (1+UiP) 
(1-u.R) 


(14a) 

(14b) 


and,  for  unity  gain  in  the  low  passband, 

■Further  details  in  the  manipulation  are  found  in  Appendix  C. 
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K  = 


(1-u. R) (1-u.R)  ...  (1-u  R) 
i  /  n 


1 

G(l) 


9 


(14c) 


where  are  the  poles  of  the  normalized  n-th  order  analog 
Butterworth  lowpass  filter  G^pjj(s)  described  in  Appendix  Al. 

Remark  2: 


Note  that  the  digital  filter  is  explicitly  dependent  on  the  factor  R. 
The  critical  frequency  will  be  determined  by  the  sampling  frequency 
through  the  ratio 

The  gain  K  may  be  arbitrarily  chosen  if  so  desired. 

A 


Example  1:  3rd  Order  Digital  Butterworth  Lowpass  Filter 

Problem:  Design  a  3rd  order  digital  Butterworth  lowpass  filter  with 
cut-off  frequency  . 

Solution:  From  Appendix  Al,  the  normalized  3rd  order  analog  Butterworth 

lowpass  filter  is  given  by 

1 

Ws)  ■  7 — — — r: — r  ■  (15a) 

(s-u^  (s-u2)  (s— u3) 

where 

n1  =  -1  (15b) 

u2  =  -. 5  +  j .866  (15c) 

u3  =  -.5  -  j .866  (15d) 
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Following  the  above  procedures: 


(a)  Select  co^/o^  =  .1476  or  tog  =  6.77^.  With  this  ratio,  the  digital 
ZOH  reproduction  of  a  sine  wave  at  is  approximately  as  shown  in 
Fig .5  . 


U) 

s 

Fig.  5 


(b)  R  =  tan(.  147 6ir)  =  .5  •  (by  choice  of  (a))  (16) 

(c)  l+(-l)(.5) 

p  _  —  .3333  (17a) 

1  1  -  (-1) (.5) 


1  +  (-.5  +  j .866) (.5) 

p  =  -  =  .  4286  +  j .  4949  (17b) 

Z  1  -  (-.5  +  j  .  866)  ( .  5) 


1  +  (-.5  -  j .866) ( . 5) 

Po  =  -  =  .4286  -  j  .4949  =  p0  (17c) 

1  -  (-.5  -  j .866) (.5)  2 
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K  = 


(.5)' 


(1  -  (-1)  ( .  5))(  1  -  (-.5  +  j.866)(.5))(l  -  (-.5  -  J  .  866)  (.5)) 


-  .04762 


(17  d) 


The  3rd  order  digital  Butterworth  lowpass  filter  is  thus  given  by 


glp(z)  ■ 


.  047  62(z  +  1)' 


(z-  .3333)  (z  -  .4286  -  j.4949)(z  -  .4286  +  j  .  4949) 


. 04762 (z3  +  3z2  +  3z  +1) 
z3  -  1.1905z2  +  .7143z  -  .1429 


(18) 


The  theoretical  frequency  response  of  G^p(z)  in  (18),  computed  as 

| G  ( e^ a)T)  |  versus  u>/ai  ,  is  shown  in  Fig.  6. 

J-*r  s 

The  recursive  equation  for  the  digital  lowpass  filter  follows  from 
(18)  as 

y(k)  =  1. 1905y(k-l)  -  .7143y(k-2)  +  .1429y(k-3) 

+  . 04762 [u(k)  +  3u(k-l)  +  3u(k-2)  +  u(k-3) ]  (19) 

where  y(k)  and  u(k)  are  respectively  the  discrete  output  and  input 
sequences  of  the  filter.  The  microprocessor-based  implementation  of  the 
lowpass  filter  given  by  (18)  or  (19)  is  described  in  Section  3.  An 
experimental  frequency  response  of  the  microprocessor-based  3rd  order 
digital  Butterworth  lowpass  filter  is  presented  in  Section  4. 
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5.1.3  DESIGN  OF  DIGITAL  HIGHPASS  FILTERS 


The  design  of  an  n-th  order  analog  Butterworth  highpass  filter 

Gjjp(s)  with  critical  cut-in  frequency  w  is  given  in  Appendix  A3.  The 

following  procedure  may  be  used  to  systematically  obtain  a  corresponding 

digital  Butterworth  highpass  filter  G  (z)  from  G  (s). 

HP 

(a)  Select  the  desired  ratio  oj,  /o>  where  ai.,  is  the  critical  cut-in 

dc  s  dc 

"r'i  v.  frequency  ;of  the  digital  highpass  f  liter  (see  Fig . 3b)  . 


_  ,  dc  , 

R  =  tan  ( -  tt) 

(0 

s 


(20a) 


so  that 


(20b) 


(c)  Using  the  conversion  scheme  (12),  a  corresponding  n-th  order 
digital  Butterworth  highpass  filter  is  obtained  as^ 


V2) 


ghp(s) 


u  - 
ac  z— 1 

R  z+1 


K(z-l)1 


(z-p1) (z— p2)  ...  (z-p  ) 


=  KG(z) 


(21a) 


with 


Details  given  in  Appendix  C. 
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1  +  u±R 


(21b) 


x  1  -  u.R 
x 

and,  for  unity  gain  in  the  high  passband, 

1 


(1  -  u  P)  (1  -  u.R)  ...  (1  -  u  R) 
1.2  n 

i:. 


G(-l)  * 


(21c) 


where  u^  are  the  poles  of  the  normalized  n-th  order  analog  Butterworth 
lowpass  filter  described  in  Appendix  Al. 

Remark  2  similarly  applies  to  the  above  design  of 

digital  Butterworth  highpass  filter. 


Example  2:  3rd  Order  Digital  Butterworth  Highpass  filter 


Problem:  Design  a  3rd  order  digital  Butterworth  highpass  filter  with 
cut-in  frequency  co^c . 

Solution:  From  Appendix  A3,  the  3rd  order  normalized  analog  Butterworth 
lowpass  filter  is  given  by 


<Ws) 


3 

s 


(s  -  u^)(s  -  u2)(s  -  Ug) 


(22a) 


where 

U;L  =  -1  (22b) 

u2  =  -.5  +  j.866  (22c) 


u^  =  -. 5  -  j .866 


(22d) 
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Following  the  procedures  outlined  above: 


(a)  Select  to  /w  =  .05  or  u  =  20  .  With  this  ratio,  the  digital 

Cl  C  S  S  ("i  c 

ZOH  reproduction  of  a  sine  wave  at  is  approximately  as  shown  in 

Fig. 7. 


(b)  P  =  tan( .  05tt)  =  .1585 


(23) 


1  +  (-1) ( .1585) 

(c)  p  =  -  =  .72637  (24a) 

1  1  -  (-1)(.  1585) 


1  +  (-.5  +  j .866) ( . 1585) 

p  =  -  =  .82366  +  j. 23195  (24b) 

1  -  (-.5  +  j .866) ( .1585) 


1  +  (-.5  -  j .866) ( . 1585)  *  . 

p„  =  - - - -  .82366  -  j. 23195  =  p„  (24c) 

J  1  -  (-.5  -  j .866) (.1585)  1 


Hence, 
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ghp(z)  ■ 


K(z  -  l)3 

(z  -  .72637)  (z  -  .82366  -  j. 23195) (z  -  .82366  +  j. 23195) 


13 

z 


K(z2  -  3z2  +  3z  -  1) 
2.3737z2  +  1.9288z  -  .5319 


(2  5a) 


with 


K  =  1/G(-1)  =  .72929 


(25b) 


The  theoretical  frequency  response  of  G  (z)  in  (25),  computed  as 
i  coT 

|G^p(e^U  ) |  versus  (o/w^  is  shown  in  Fig. 8. 

The  recursive  equation  for  the  digital  highpass  filter  follows  from 
(25)  as 

y (k)  «  2 . 3737y(k-l)  -  1.9288y(k-2)  +  ,5319y(k-3) 

+  .72929 [u(k)  -  3u(k-l)  +  3u(k-2)  -  u(k-3)]  ,  (26) 

where  y(k)  and  u(k)  are  respectively  the  output  and  input  sequences  of  the 
filter.  The  microprocessor-based  implementation  and  experimental 
frequency  response  of  the  highpass  filter  (25)  or  (26)  are  described  in 
Sections  3  and  4. 
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5.1.4  DESIGN  OF  DIGITAL  BANDPASS  FILTERS 


The  design  of  a  2n-th  order  analog  Butterworth  bandpass  filter 
G  (s)  with  bandwidth  BW  and  mid-band  frequency  to  is  given  in  Appendix  A5, 

x5  r  3.0 

From  G^pCs) ,  a  corresponding  digital  Butterworth  bandpass  filter  ^-gpC2) 
can  be  obtained  as  follows: 


(a)  Choose  to,  /to  ,  where  to,  is  the  mid-band  frequency  of  the  digital  filter, 

do  s  do 

so  that  (to,  +  —x~)  /to  <<  .5. 
do  2  s 

(b)  Obtain 


to, 

„  ,  do  . 

R  =  tan  ( -  7 r) 

to 

s 


to  T 
ao 


2 


or 


2F 

to 

ao 


(27a) 


(27b) 


(c)  Invoking  the  conversion  scheme,  a  corresponding  2n-'th  order  digital 

1  2 

Butterworth  bandpass  filter  can  be  obtained  as  ’ 


Vz)  ■  Vs) ' 

:  2  z— i 

I  S  T  Z+l 


K(z  -  l)n(z  +  l)n 


*  * 
(z  -  p.)(z  -  Pl)  ...  (z  -  p  )(z  -  p  ) 
li  n  n 


(28a) 


*  denotes  complex  conjugation. 

Details  of  manipulation  is  given  in  Appendix  C. 
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with 


P 


i 


A 


(2/t  +  cj 

(2/t  -  c .) 

1 


i~  1  j  •  •  * )  H) 


(28b) 


BW 

2 


u. 

l 


+  JO) 


ao 


(28c) 


(BW)n(2/r)n 
(2/r  -  c.)(2/t  -  c.) 

l  l 


(28d) 


where  t  is  given  by  (27). 


Remark  3 :  The  design  of  the  bandpass  filter  (28)  assumes  that 

2  2 

|BWu_L|  «  (see  Appendix  A5)  .  This  is  equivalent  to  considering  a 

bandpass  filter  with  a  high  Q-f actor  (the  ratio  of  the  midband  frequency 
to  the  bandwidth),  i.e.. 


*  m 

0  =  -^r  >  1. 

BW  — 


(29) 


A 
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Example  3:  6th'  Order  Digital  Butterworth  Bandpass  Filter 


Problem:  Design  a  6th  order  digital  Butterworth  bandpass  filter  with 

bandwidth  of  BW  =  2tt  x  20  rad/s  and  midband  frequency  of 

to,  =  2x  x  20  rad/s.  (Note  the  Q-factor  =  =  1.) 

°  BW 

Solution:  From  Appendix  A5,  a  6th  order  analog  Butterworth  bandpass 


filter 

having  the  above  specification  (BW  =  2tt  x  20, 

to  = 

ao 

io^o  =  2-rrx  20)  is  given  by 

3  3 

(BW)  s 

gbp(s)  = 

&  it  * 

(s  -  c1)(s  -  c1)(s  -  c2)(s  -  c2)(s  -  C3HS  -  c^) 

(30a) 

with 

* 

(2tt  X  20)  (-1) 

cr  °i  ■ 

-  +  j2ir  x  20 

2 

= 

-20tr  +  jAOir 

(30b) 

(2tt  x  20) 

C2’  °3 

=  -  (- .  5  +  j  .  866)  +  j  2tt x  20 

2 

=  -10tt  +  j57.32n,  -10tt  -  j22.68TT 

(30c) 

*  * 
c2  *  c3 

=  -10tt  -  j57.32ir,  -10tt  +  j22.68iT 

(30d) 

(a)  Select  to,  /to  =  .2  or  to  =  5<o,  .  With  this  ratio,  the  relative 
do  s  s  do 

position  of  the  passband  with  respect  to  the  sampling  frequency  is 
approximately  as  shown  in  Fig.  9. 
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H  |tO 


sampling 

frequency 


Fig.  9.  BandPass  Filter 


(b)  R  =  tan(.2rr)  =  .73 


=  1 17.2 


*  172  +  (-20tt  +  j40tt)  4? 

(c)  p.f  p  =  -  =  .6246e^jl'J4/b 

X  1  172  -  (-20tt  +  j 40ir) 


=  .1382  +  j . 6091  «  .14  +  j . 61 


*  172  +  (-107r  ±  j57.32ir)  _ +11.6326 

,  p0  =  -  =  .  8409 e—J 


2  ’  v2 


172  -  (-IOtt  +  j 57 . 32ir) 


=  -.0519  +  j .8393  «  -.052  +  j.84 


*  172  +  (-10,  +  j  22. 68,)  _+1.8194 

p„,  p-  =  -  =  ,7319e-J 

J  172  -  (-IOtt  +  j22.68^) 


=  .4996  +  j . 5348  =  .50  +  j  .  53 


(32c) 


2  2  *  2 

(x  "  C1)(T  "  cl)  =  266-35 

2  2  *  2 

(T  '  C2)(7  "  c2}  =  271*68 

2  2  *  2 

(f  -  c_)(-  -  cj  =  215.54^ 

T  3  T  3 

(2tt  x  20)3(172)3 

K  =  - — - — - —  =  .0415  (32d) 

266.35  x  271.68  x  215.54 

Hence, 

0.04l5(z  -  1) 3  (z  +  l)3 

CL„(z)  -  - - 

TJp 

(z  -  .14  -  j . 61) (z  -  .14  +  j . 61) (z  +  .052  -  j.84)(z  +  .052  +  j .84) 
1 

x  - 

(z  -  . 5  -  j . 53) (z  -  .5  +  j.53) 

0.0415(z6  -  3z4  +  3z2  -  1) 

(z6  -  1.176z5  +  1.7778z4  -  1.3219z3  +  1.0035z2  -  .3611z  +  .1473) 

(33) 

The  theoretical  frequency  response  of  G^^(z)  in  (33) ,  computed  as 
|Ggp(e3a)rr)  |  versus  cj/cos,  is  shown  in  Fig.  10. 

The  recursive  equation  for  the  6th  order  digital  bandpass  filter  follows 
from  (33)  as 

y(k)  =  1.17 6y(k-l)  -  1.7778y(k-2)  +  1.3219y(k-3)  -  1.0035y(k-4) 
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+  .3611yCk>5)  -  .  147 3y  (k-6)  +  ..04l5u(k)  -  ..1245u(k-2)  +  ,1245u(k-4) 
-  . 0415u(k-6)  ('34) 


where  y(k)  and  u(k)  are  respectively  the  output  and  input  of  the 

filter.  The  microprocessor-based  implementation  and  the  experimental 
frequency  response  of  the  bandpass  filter  (33)  or  (34)  are  given  in  Sections 
3  and  4. 
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5.1.5  DESIGN  OF  DIGITAL  BANDSTOP  FILTERS 


The  design  of  a  2n-th  order  Butterworth  bandstop  filter  G^-Cs)  with 

bandwidth  BW  and  midband  frequency  to  is  given  in  Appendix  A6.  Using 

G  (s)  ,  a  corresponding  digital  Butterworth  bandstop  filter  G  (z)  can  be 
Bo  BP 

obtained  as  follows: 

(a)  Choose  to,  / to  ,  where  to,  is  the  desired  midband  frequency  of  the 

do  s  do 

digital  filter,  so  that  (to,  +  ~x-)/to  <<  .5. 

do  2  s 

(b)  Obtain 


R  = 


to, 

tan  ( — —  tt) 

“s 


(0  T 

ao 


(35a) 


so  that 


2R 

to 

ao 


(35b) 


(c)  Using  the  bilinear  transformation,  a  corresponding 2n-th  order 
digital  Butterworth  bandstop  filter  can  be  obtained  as 


gbs(z) 


gbs<s>  i 


S  = 


2  Z— 1 
t  z+1 


K(z  -  z  )n(z  -  z  ^)n 
o  o 


(z  -  p^(z  -  p^  ...  (z  -  Pn)(z  -  Pn) 


*  » 


(36a) 


where 


pi " 


4  2/T  +  ci 


2/t  -  c± 


(36b) 


A  BW 

C  .  =  -r-  U.  +  JtO  , 

i  2  i  —  J  ao 


(see  Appendix  A5) 


(36c) 
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K  = 


[  (2/t)  2  +  a)2  ]n 
ao 


(2/t  -  c- ) (2/t  -  c  )  ...  (2/t  -  c  ) (2/t 
11  n 


* 

Cn) 


Zcr  =  (2/t  -  jo)ao)/(2/T  +  j(oaQ)  . 


Example  4:  6th  Order  Digital'  Butterworth  Bandstop  Filter 


(36e) 

(36f) 


Problem:  Design  a  6th  order  digital  Butterworth  bandstop  filter  with 

BW  =  2ir  x  20  rad/s  and  midband  frequency  of  =  2tt  x  20  rad/s. 

(Note  the  0-factor  =  u,  /BW  =  1.) 


Solution:  From  Appendix  A6,  a  6th  order  analog  Butterworth  bandstop 

filter  having  the  above  specification  (BW  =  2ir  x  20  and 

a)  =  to,  =  2tt  x  20)  is  given  by 
ao  do  1 

r  2  ±  2  ^ 

(s  +  0)  ) 

G  (s)  =  - - - §2 - - - —  (37a) 

(s  -  c1)(s  -  Cj)(s  -  c2)(s  -  c2)(s  -  c3)(s  -  c3) 


cr 

c^  =  -20tt  +  j40ir 

(37b) 

* 

C2  * 

c2  =  -10ir  +  j  57 . 32tt 

.  (37c) 

* 

C3’ 

c3  =  -10tt  +  j22.68r 

(37d) 

We  note  that  c^  in  (37)  are  the  same  as  those  of  (30)  due  to  certain 
similarities  in  the  design  specification  of  Examples  3  and  4.  For 
convenience,  let  us  also  choose  the  design  variables  as  in  Example  3. 
That  is: 

(a)  tojo/<i>s  =  .2.  With  this  ratio,  the  relative  position  of  the  stopband 
with  respect  to  the  sampling  frequency  is  approximately  as  shown  in 
Fig.  11  , 
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(f-  C^cf  -  c*)  =  266. 352  , 

(f  “  C2)(f  "  c2)  =  271-682  *  (40) 

(f  -  c  )(-  -  c*)  =  215. 542  . 

T  J  T  3 


In  addition,  we  compute 


(-)  +  a)2  =  1722  +  (2tt  x20)  2  =  45375 

t  ao 

so  that  (36e)  yields 


(41a) 


K  = 


45375' 


266. 352  x  271. 682  x  215. 542 


=  .384 


(41b) 


The  zeroes  specified  by  (36d)  are  given  by 


172  +  j2ir  x  20  _  ^  2619 

z  =  -  =  e  J  =  , 3040  _  j  .9527  (42a) 

°  172  -  j 2-rr  x  20 


z'1  =  .3040  +  j .9527  . 

o  J 


(42b) 


From  the  above,  the  corresponding  6th  order  digital  Butterworth  bandstop 
filter  is  obtained  as 


gbs(z)  ■ 


,384(z  -  .3040  +  j.9527)3(z  -  .3040  -  j.9527)3 
(z  -  .14  -  j . 61) (z  -  .14  +  j . 61) (z  +  .052  -  j.84)(z  +  .052  +  ; 


(z  -  .5  -  j  .53)  (z  -  .5  +  j.53) 


.84) 
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(43) 


. 384 (z6  -  1.824z5  +  4.109z4  -  3.873z3  +  4.109z2  -  1.824z  +  1) 
z6  -  1.176z5  +  1.7778z4  -  1.3219z3  +  1.0035z2  -  .3611z  +  .1473 


The  theoretical  frequency  response  of  G  (z)  in  (43),  computed  as 

|  *  ^  t  "  ^ 

iG  (eJ  )|  versus  u/u)  is  shown  in  Fig.  12. 

I  ij  o  •  S 

The  recursive  equation  for  the  6th  order  Butterworth  bandstop  filter 
follows  from  (43)  as 

y (k)  =  1. 176y(k-l)  -  1. 7778y(k-2) •+  1.3219y(k-3)  -  1.0035y(k-4) 

+  . 3611y(k-5)  -  . 1473y(k-6)  +  .384[u(k)  -  1.824u(k-l)  +  4.109u(k-2) 

-  3 . 873u(k-3)  +  4 ,109u(k-4)  -  1.824u(k-5)  +  u(k-6)],  (44) 


where  y(k)  and  u(k)  are  respectively  the  output  and  input  of  the  digital 
filter.  The  microprocessor -based  implementation  and  the  experimental 
frequency  response  of  the  bandstop  filter  (43)  or  (44)  are  given  in 
Sections  3  and  4. 
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5] 


5.1.6  A  FURTHER  EXAMPLE 


The  effect  of  pole  variation  in  the  filter  may  be  vividly  illustrated 
by  considering  the  following  example. 

Example  5:  3rd  Order  Digital  Chebychev  Lowpass  Filter 

Problem:  Design  a  3rd  order  digital  Chebychev  lowpass  filter  with  cut-off 
frequency  u(jc»  having  a  ripple  factor  r  =  1  db.  Similar  design 
considerations  for  Example  1  may  be  used. 

Solution :  Let  a  3rd  order  analog  Chebychev  lowpass  filter  be  denoted  by 

K 

G  (s)  =  - § -  (45) 

LPC  (s  -  u’)(s  -  u’)(s  -  u’) 

where  the  poles  u^,  and  u^  may  be  determined  as  follows. 

The  specification  of  the  filter  translates  into  the  frequency  response 
sketched  in  Fig.  13..  The  ripple  factor 

r  =  1 

1 

=  20  log(l)  -  20  log(  t - 5-  ) 

/  1  +  e 

=  10  log (  1  +  e2  ),  (46a) 

so  that 

e  = 

=  .5088  , (r  =  1)  .  (46b) 
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Fig.  13.  Chebychev  Lowpass  Filter 


Let  n  =  3  denote  the  order  of  the  filter  and  define 


a 


—  sinh 
n 


-  In  (-+/(-) 2  +  1  ) 
n  e  v  e 


3  ln  (  .5088  V (  .5088  5  +  1  ^ 


=  .4760  , 


so  that 


tanh  a  = 


a  -a 
e  -  e 


a  -a 
e  +  e 


(47a) 


=  .4430  .  (47b) 

The  poles  of  the  analog  Chebychev  lowpass  filter  may  be  treated  as 
being  the  poles  of  the  analog  Butterworth  lowpass  filter  whose  real  parts 
are  reduced  by  a  factor  of  tanh  a.  Using  the  values  of  u.^  in  Example  1, 
the  poles  of  the  "normalized"  analog  Cheybechev  filter  (  45  )  can  hence  be 
determined  from 
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u'  =  Re  {u.}  x  tanh  a  +  ilm{u.}. 
i  1  J  1 


(49a) 


=  -.2215  +  j.866 


(49b) 


ul  =  -.2215  -  j .866 
3  J 


(49c) 


Using  relationship  (14b) ,  the  poles  u|  in  the  s-domain  can  be  mapped  into 
the  poles  in  the  z-domain  as 


Pi  = 


(1  +  u|R) 
(1  -  u IR) 


With  the  same  choice  of  R  as  in  Example  1  (i.e.,  R  =  .5),  we  obtain 


1  +  (-.443) ( . 5) 
1  _  (-.443) ( . 5) 


=  .644 


(51a) 


1  +  (-.2215  +  j . 866) ( . 5) 
1  -  (-.2215  +  j . 866) ( . 5) 


,8269ej'8247  =  .563  +j. 609  (51b) 


^  »  _  _  1 
p3  P2 


=  .563  -  j . 609 


(51c) 


Hence,  a  corresponding  3rd  order  digital  Chebychev  lowpass  filter  can  be 
obtained  as 


K(z  +  1) 


glpc(z)  ■ 


(z  -  p.p  (z  -  pp  (z  -  pp 


54 


K(z  +  l)3 

(z  -  . 644) (z  -  .563  -  j.609)(z  -  .563  +  j . 609) 

K(z3  +  3z2  +  3z  +  1) 

(z3  -  1.77z2  +  1. 413z  -  .443) 


=  KG(z) 


(52a) 


where  for  unity  gain  at  the  low  frequency 


1  .2 

K  =  -  =  —  =  .025  .  (52b) 

6(1)  8 


as 


G 


The  theoretical  frequency  response  of  G 
T  -(e^7)  |  versus  to/w  is  shown  in  Fig. 


1PC(z)  given  by  (52),  computed 
14. 


The  recursive  equation  for  the  3rd  order  Chebychev  lowpass  filter 
follows  from  (52)  as 


y (k)  =  1. 77y(k-l)  -  1.413y(k-2)  +  .443y(k-3) 

+  . 025u(k)  +  . 075u(k-l)  +  .075u(k-2)  +  .025u(k-3)  ,  (53) 


where  y(k)  and  u(k)  are  the  output  and  the  input  of  the  digital  filter. 

The  microprocessor-based  implementation  and  experimental  frequency 
response  of  the  digital  Chebychev  lowpass  filter  (52)  or  (53)  are  given  in 
Sections  3  and  4. 
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5.2  MICROPROCESSOR  REALIZATION  OF  DIGITAL  FILTERS 


The  microprocessor  firmware  for  implementing  the  digital  frequency 
selective  filters  designed  in  Section  2  is  described  below. 

Fig. 15  shows  the  block  diagrams  of  the  hardware  for  the  microprocessor- 
based  signal  processing  system  used  in  the  realization  of  the  digital  filters. 
The  CPU  of  the  system  is  the  Motorola  MC6802  microprocessor,  which  is 
supported  by  2716  EPROM  (monitor  program),  6810  RAM's,  6821  PIA's,  keypad, 
L.E.D.  display  and  interfacing  buffers  to  form  the  microcomputer  called 
MOUSE'*'  (see  Fig.  16  ).  The  microcomputer  operates  at  a  clock  rate  of  1  MHz. 

The  data  acquisition  unit  is  the  Datel  MDAS-16  multiple  (multiplex) 
channel  12-bit  A/D  converter  with  a  conversion  time  of  about  20  ys  per  data. 

A  Datel  Hzl2BGC  12  bit  D/A  converter  with  a  settling  time  of  3  ys  is  used 
as  a  zero-order-hold  output  of  the  microprocessor-based  system.  The  inter¬ 
face  between  the  8-bit  MC6802  and  the  12-bit  I/O  (A/D  and  D/A)  peripherals 
are  done  through  6821  PIA's  (Fig. 17).  The  software  (<  lk  bytes)  for  the 
digital  filters  are  stored  in  the  external  2716  EPROM.  The  additional 
2114  RAM's  provide  handy  facilities  for  debugging  and  immediate  alterations 
of  the  software  if  desired.  The  wiring  diagrams  for  the  microprocessor- 
based  signal  processing  system  is  shown  in  Figs.  16  and  17. 


Acronymn  for  Microcomputer  of  Oakland  University  School  of  Engineering. 
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The  main  consideration  in  the  microorocessor  software  for  the  digital 
filters  involves  the  development  of  fast  and  efficient  arithmetic  subroutines, 
the  scaling  of  the  recursive  digital  filter  equations,  the  handling  of 
saturation  and  some  memory  management. 

A  fast  3  bytes  x  1  byte  multiplication  subroutine  with  an  execution 
time  of  about  96  ps  was  developed  for  the  digital  filter  implementation. 

Other  main  subroutines  include  a  3  bytes  +  3  bytes  summation 

(about  85  ps) ,  transfer  and  negation  of  3  byte  data.  Details  of  these 
subroutines  are  given  in  Appendix  B. 

In  order  to  minimize  the  occurrence  of  saturation  (overflow  or  underflow) 
in  the  finite  wordlength  data  (3  bytes  or  24  bits),  the  recursive  formulas 
for  each  of  the  filters  will  be  scaled  in  a  fashion  similar  to  those  done 
in  an  analog  computer  simulation.  The  scaled  recursive  equations  for  the 
digital  filters  from  Examples  1-5,  are  shown  in  Table  1. 

It  is  remarked  that  the  forward  gains  K  of  the  filters  may  be  reduced, 
if  necessary,  to  achieve  a  proper  scaling  which  will  not  saturate  the  filter 
output.  Alternatively,  the  output  saturation  can  be  handled  by  use  of 
overflow  test  instructions  in  the  microprocessor  software. 

The  6802  microprocessor  software  for  implementing  the  digital 
frequency  selective  filters  described  in  Table  1  are  given  in  Appendix  B. 
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r 


Fig.  15  Block  Diagram  of  Microprocessor  System 
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Fig. i7  Microprocessor  Controlled  ADC  4  DAC 


Table  1:  SCALED  RECURSIVE  EQUATIONS  FOR 


THE 


DIGITAL  FILTERS 


Butterworth  Low-Pass  (Example  1) 

. 84y(k)  =  y(k-l)  -  .6y(k-2)  +  .12y(k-3) 

+  . 04u(k)  +  . 12u(k-l)  +  . 12u(k-2)  +  .04u(k-3) 


Butterworth  High-Pass  (Example  2) 

. 4213y(k)  -  y(k-l)  -  .8l26y(k-2)  +  .2241y(k-3) 

+  .3072u(k)  -  . 9217u(k-l)  +  .9217u(k-2)  -  .3072u(k-3) 


Butterworth  Band-Pass  (Example  3) 

. 5625y(k)  =  . 6615y(k-l)  -  y(k-2)  +  .7436y(k-3)  -  .5645y(k-4)  +  .2031y(k-5) 
-  .0829y(k-6)  +  .0233u(k)  -  .0700u(k-2)  +  .0700u(k-4) 

. 0233u(k-6) 


Butterworth  Band-Stop  (Example  4) 

. 5625y(k)  =  ,6615y(k-l)  -  y(k-2)  +  .7436y(k-3)  -  .5645y(k-4)  +  .2031y(k-5) 

-  . 0829y(k-6)  +  .216u(k)  -  .394u(k-l)  +  .8875u(k-2) 

-  . 8366u(k-3)  +  . 8875u(k-4)  -  .394u(k-5)  +  .216u(k-6) 


Chebychev  Low-Pass  (Example  5) 

. 565y(k)  =  y(k-l)  -  .7983y(k-2)  +  .2503y(k-3) 

+  . 0l41u(k)  +  . 0423u(k-l)  +  .0423u(k-2)  +  .0l41u(k-3) 
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5*  3  EXPERIMENTAL  FREQUENCY  RESPONSE  OF  MICROPROCESSOR-BASED  DIGITAL  FILTERS 


The  experimental  set-up  for  recording  the  frequency  response  of  the 
microprocessor-based  digital  filters  is  shown  in  Fig.  18  .  The  test  input 
u(t),  generated  by  the  voltage  controlled  oscillator,  consists  of  a  constant 
amplitude  sinusoidal  signal  whose  frequency  modulates  or  sweeps  (sufficiently) 
slowly  from  low  frequency  to  high  frequency  and  vice  versa,  i.e. , 

u(t)  =  A  sin  ait 

where  the  frequency  a)  is  controlled  by  a  triangular  or  saw-tooth  signal  w(t) . 
The  digital  output  of  the  filter  is  recorded  on  a  storage  scope  whose 
horizontal  axis  is  driven  by  the  same  w(t).  From  the  set-up,  one  can 
experimentally  determine  the  frequency  responses  of  the  microprocessor-based 
digital  filters.  Figs.  19-23  show  the  actual  experimental  frequency 
responses  of  the  microprocessor-based  3rd  order  filters  designed  in  the 
examples  of  Section  2. 

For  comparison,  the  critical  frequencies  of  the  theoretical 
filters  and  the  implemented  microprocessor-based  filters  are  tabulated  in 
Table  2.  As  shown  in  the  table,  the  specification  of  the  filters  in 
terms  of  and  BW  have  been  met  satisfactorily. 

It  is  important  to  note  that  the  critical  frequencies  and 

can  readily  be  altered  by  simply  adjusting  the  sampling  frequency  id  . 

s 

There  is,  however,  an  upper  bound  on  the  maximum  possible  sampling  frequency 

which  can  be  used  for  the  filter  due  to  the  finite  speed  of  the  microprocessor. 

Nevertheless,  the  design  specifications  concerned  in 

the  present  investigation  can  be  satisfactorily  fulfilled  by  the  current 

generation  of  8-bit  microprocessors.  For  more  stringent  design  specifications, 

one  may  resort  to  the  new  generation  of  16-bit  microprocessors  and/or 

use  of  high  speed  arithmetic  logic  chips. 


u(t)  =  A  sin(a>(t)+«0 
w(t)  ~W(t) 

Voltage 
Controlled 
Oscillator 

(Digital  Filters) 


Microprocessor 

-Based 

Signal  Processor 


Fig. 19  Experimental  Frequency  Response  of  Microprocessor- 

Based  3rd  Order  Butterworth  Lowpass  Filter  Gjp(z)>  (Example  1) . 


Fig. 20  Experimental  Frequency  Response  of  Microprocessor-Based 

3rd  Order  Butterworth  Highpass  Filter  G  (z)  (Example  2) 

HP 
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Fig *21  Experimental  Frequency  Response  of  Microprocessor— Based 

6th  Order  Butterowrth  Bandpass  Filter  G  (z )  (Example  3) 

BP 


Fig. 22  Experimental  Frequency  Response  of  Microprocessor-Based 
6th  Order  Butterworth  Bandstop  Filter  G  (z)  (Example  4) 

DO 
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Fig. 23  Experimental  Frequency  Response  of  Microprocessor^Based 
3rd  Order  Chebychev  Lowpass  Filter  (Example  5) 


5. 
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Table  2.  Critical  Frequencies  of  Theoretical  and  Microprocessor-Based  Filters 


Digital 

Filter 

Sampling 

Frequency 

03  (Hz) 
s 

Cut-Off  or 

Cut-In 

Frequency 

“dc  (Hz) 

Frequency  at 

Attenuation 

=  0.1 

“o.i  (Hz) 

Remarks 

3rd  Order 

Butterworth 

Lowpass 

Theoretical 

200 

30 

52 

Experimental 

200 

27 

49 

3rd  Order 

Chebychev 

Lowpass 

Theoretical 

200 

29 

44 

The  cut-off  rate  is 

faster  than  that  of 

the  Butterworth  filter. 

Experimental 

200 

27 

45 

3rd  Order 

Butterworth 

Highpass 

Theoretical 

200 

10 

5 

Experimental 

200 

9.5 

5 

Digital 

Filter 

Sampling 

Frequency 

o)  (Hz) 
s 

Midband 

Frequency 

oj  (Hz) 
o 

First 

Aliasing 

Midband  Freq. 
03^  (Hz) 

■ 

Bandwidth 

BW  (Hz) 

Remarks 

6th  Order 

Butterworth 

Bandpass 

Theoretical 

100 

25 

75 

mm 

03/03  =  .25 
o  s 

0)  /  U3  =  .2 
o  s 

Experimental 

140 

28 

100 

17 

6th  Order 

Butterworth 

Bandstop 

Theoretical 

100 

20 

80 

14 

- 

Experimental 

104 

16 

78 

16 
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APPENDIX  A 


SUMMARY  OF  ANALOG  BUTTEPWORTH  FREQUENCY  SELECTIVE  FILTERS 


Al.  Low-Pass  Filter  G^p^Cs)  with  Normalized  Cut-Off  Frequency ' (“ac  =  1) : 


glpn(s)  " 


(s  -  u1) (s  -  u2> 


(s  -  un) 


where  u^,  i  =  1,  . . . ,  n,  are  the  stable  poles  which  lie  on  the  unit  circle 
in  the  s-plane  as  shown  in  Fig.Al. 


n  =  odd  (e.g.,  n=3)  n  =  even  (e.g.,  n=4) 


Fig.  Al.  Pole  Locations  of  Normalized  Butterworth  Lowpass  Filter 

^  We  note  that  complex  poles  must  occur  in  complex  conjugate  pairs  in  order 
for  the  filter  to  be  physically  realizable. 


A.  1 


A2:  Low-Pass  Filter  G^(s)  with  Arbitrary  Cut-Off  Frequency  o)  : 


To  translate  GT_„(s)  into  G.^Cs),  one  substitutes  s/03  for  s  in 
LPN  LP  ac 

(A. 1)  and  obtain 


03 


GLp(s)  = 


ac 


03 


ac 


(s  -  “acV  (s  '  “ac11!5  ** 


A  typical  set  of  pole  locations  of  G. 


LP 


.  (s  -  03  u  )  (A. 2) 

ac  n 

(s)  is  shown  in  Fig.A2. 


Fig.  A2.  Pole  Locations  of  Butterworth  Lowpass  Filter  with  03  (n  =  5) 

ac 


A3:  High-Pass  Filter  G^^Cs)  with  Normalized  Cut-In  Frequency  ( to  c  =  1)  : 


To  translate  G^p^(s)  into  G^p^(s),  one  substitutes  1/s  for  s  in 
(A.  1)  and  obtains 

s  s  s 

^npM  ^  S  ^  ~  ... 

(1 — su  )(1  -  su  )  (1  -  su  ) 

12  n 

(-l/u1)s  (-l/u2)s  (-l/un)s 

.  •  •  •  ——————— 

(s  -  1/u.)  (s  -  1/u.)  (s  -  1/u  ) 

12  n 

s  s  s 

_  ^ ...  ,  (A. 3) 

(s  -  u.,)  (s  -  u_)  (s  -  u  ) 

1  2  n 

where  the  last  equality  follows  from  the  fact  that  1/u^  =  u^,  |u^|  =1 

and  u^  occur  in  complex  conjugate  pairs.  A  typical  set  of  poles  and  zeroes 
for  Gjjpjy ( s)  is  shown  in  Fig. A3. 


Fig.  A3. 


Pole-Zero  locations  of  Normalized  Butterworth  Highpass  Filter,  (n 


A4:  High-Pass  Filter  G  (s)  with  Arbitrary  Cut-In  Frequency  <u  : 

Hr  B.C 


To  translate  G,TT,..(s)  to  GT(s)  ,  one  substitutes  s/w  for  s  in  (A. 3) 
HPN  HP  ac 


and  obtains 


°HP(S)  = 


(s  -  0)  u.)  (s  -  0)  U-)  (s  -  to  u  )  , 

ac  1  ac  2  ac  n 


where  a  typical  set  of  poles  and  zeroes  for  G  (s)  is  shown  in  Fig.A4. 

Hr 


Fig.  A4.  Pole-Zero  Locations  of  Buttetworth  Highpass  Filter  with  ui  .  (n  =  4) 

ac 


A.4 


A5:  Band-Pass  Filter  G^Cs)  with  Bandwidth  BW  and  Midhand  Frequency  oi  : 

BP  ^  J  ao 


A  high  Q-factor  will  generally  be  assumed,  i.e., 
To  translate  G^p^(s)  to  G^Cs),  one  substitutes 


0  =  0)  /BW 
ao 


1. 


1 

BW 


2 

s 


+  0) 


ao 


for  s  in  (A.l)  and  obtains 


BW.  s 


BW.  s 


BW.  s 


rbp(s)  =  ;  2 


(s" — ;BWu..  s  +  0)2  )  (s2  -  BWu.s  +  a)2  )  (s2  -  BWu  s  +  to2  ) 

-  _  1  ao  2  ao  n  ao 


BW.  s 


BW.  s 


BW.s 


(s  -  PjXs  -  qx)  (s  -  p2)(s  -  q2>  (s  -  Pn>(s  -  q_) 


BW.s 


BW.s 


BW.s 


^  ^  ^ 
(s  -  c1 )  Cs  -  C.)  (s  -  C— )  ( s  -  cj  (s  -  c  )  ( s  -  C  ) 

11  2  2  n  n 


,  (A. 5) 


where 


>i>  q±  =  ^  ±  %/( BWuJ2  - 


4oj 


ao 


which  may  be  approximated  by 


*  A  BW 

V  ci  =  T  ui  ±  3“ao 


for  a  high  0-factor.  A  typical  pole-zero  location  of  G  (s)  is  shown  in 

or 

Fig.  A5. 


A.  5 


m/2 


Fig.  A5. 


Pole-Zero  Locations  of  Butterworth  Bandpass 
Bandwidth  BW  and  Midband  Frequency  u)  q. 


Filtei 


Fig.  A6.  Pole-Zero  Locations  of  Butterworth  Bandstop  Filtei 

n  t  t»t t j  t/j  1  T-i 


[§|i] 


A6:  Band-Stop  Filter  Ggg(s)  with  Bandwidth  BW  and  Midband  Frequency  co  : 


A  high  O-factor  will  similarly  be  assumed,  i.e.,  0  =  co  /B W  >  1. 

ao  — 

To  translate  G,T_,T(s)  to  G^Cs),  one  substitutes 
HPN  BS 


BW 


2  -  2 
S  +0) 

ao 


for  s  in  (A. 3)  and  obtains 


rJBS(s)  ,  2 


,  2  .  2  . 

(S  +  CO  ) 

ao 


,  2  ^  2  \ 
(s  +  CO  ) 
ao 


(s^  -  BWu-S  +  co2  )  (s2  -  BWu„s  +  co2  ) 
l  ao  l  ao 


,2  2> 
(s  +  CO  ) 
ao 


(s2  -  BWu  s  +  co2  ) 
n  ao 


(s+jcoao)  (s-jcoao)  (s+jcoao)  (s-jcoao) 


*  *  *  * 
(s  -  c^Cs  -  c^)  (s  -  c2)(s  -  c2) 


(s+jcoao)(s-jcoao) 

* 

(s  -  c  ) (s  -  C  ) 
n  n 


(A.  6) 


where  c^  are  as  defined  for  (A5) .  A  typical  pole-zero  location  for  Ggg(s) 
is  shown  in  Fig.  A6. 
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MICROPROCESSOR__PROGRANS_FOR  DIGITAL  FREQUENCY  SELECTIVE  FILTERS 
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APPENDIX  C 


DERIVATION  OF  DIGITAL  FILTERS 


C.l.  Derivation  of  Eq.  (14) 
From  (A. 2), 


p(s)  - 


dC  •  •  • 


[S— (i)  ll1  )  (s-o)  U V)  (  S— a)  u  ) 

v  ar  I  '  '  o r*  J*  '  a/-*  1 


Using  the  substitution  in  (14),  one  obtains 


glpU) 


GLP(s)  a) 


ac  z-1 
R  z+I 


ac  z-1 
R  z+I 


/  ac  z-1  x 

(1 r  z+T  ■  “acun  } 


R(z+1)  ...  R(z+1) 

.z-1  -  Ru1(z+l)j  [z-1  -  Run(z+1). 


u+i: 

1+Ui  R 


l+u  R 


(1  ■u1R)...(l-unR)  (z  i_u|R)***(z  -  l_unR^ 


U+D 

;z-p1)...(z-pn) 


(c.f.  Section  2.2) 


C.2.  Derivation  of  Eq.  (21) 
From  (A. 3), 


GHP(S)  "  (s-»acUl)  (s-o,acu2)  (s-^acun) 


C.l 


Using  the  substitution  in  (21),  one  obtains 


Gup(z)  -  Gud (s) 


u  n 

ac  z-1 
R  z+T 


U)  >1 

ac  z-1 
R  z+1 

/“ac  z-1  _  v 
(  R  z+1  “acV 


UJ  1 

ac  z-1 
R  z+1 

Q)  i 

/  ac  z-1  > 

V — TT - XT"  “  <i>  U  ) 

R  z+1  ac  n ' 


z-1 _ 

[z-1  -  u^R(z+l); 


_ z- 1 _ 

[z-1  -  unR(z+l)] 


(l-u1R)...(l-unR)  (z  - 


(z-1) 

i+UjR 

I^TR‘}  (z 


=  K  Tr-rr^T— ^ 


k z  PiJ***(z-Pn- 


(c.f.  Section  2.3). 


C.3.  Derivation  of  Eq.  (28) 
From  (A. 5) 


r  /_\  _  BW.s _  BW.s _ 

BP  (s-c-i ) (s— Ci *)  (s-cl)(s-c,  *) 


Using  the  substitution  in  (28),  one  obtains 

GBP (2)  =  GBP^  2  z  1 
s  =  t  'FT 


gw  —  --A 
BW  t  z+1 


(1  Zzi_  c  )(2H  _  c  *) 

't  Z+1  clMt  z+1  C1  ' 


gW  —  iLlL 

X  Z+1 _ 

/  2  z— 1  w  2  z— 1 

(  T^T  "  V(TFT 


BW  I-  ( z— 1 )  ( z+1 ) 


(rci>  (rci  > 


tfV 


BH  |  (z-1) (z+l) 


iAi  ,2  ,  h2,  *w.  <t**nJ,r.  (tcn*> 


(rc„»-K 


C.2 


A  K  (Z-I)n(z+I)n 

( z-px )  ( z-p^) . . .  ( z-pn )  ( z-pn* ) 


(c.f .  Section  2.4) . 


C.4.  Derivation  of  Eq.  (36) 
From  (A. 6), 


(stj»aoKs-j«.ao) 


BS 


<»>  •  (s-c^fTs-c^r 


(s-j^0)(s-j.a0) 

(s-cn)‘s'cn*) 


Using  the  substitution  in  (36),  we  obtain 


GBS^  _  GBS^ 


2  z— 1 
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z+1  J<ilao^t  z+1  “  Jwao^ 


-2  z-l 


,2  z-l 
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T  ~J“ao 


z  +^Mao 
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